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We  discuss  the  numerical  solution  of  the  nonlinear  Volterra  integrodif ferential 
equation 


t 

yx'(t)  + J a(t-s)F(x(t)  , x(s))ds  = 0 , t >_  0 , 


x(s)  = g(s)  , “«>  < s < 0 


Here  x(t)  is  the  unknown  function,  g is  given  history,  y is  a small  positive 
number,  a(t)  is  a positive,  decreasing,  logarithmically  convex  kernel,  while  F 
vanishes  on  the  diagonal,  and  is  increasing  in  the  first  and  decreasing  in  the 
second  variable.  The  time  discretization  is  done  using  backwards  differences,  and 
we  show  that  the  discretization  preserves  the  qualitative  properties  of  the  solu- 
tions and  give  an  error  bound  which  is  uniform  in  y and  t for  0 < y <_  y^, 
and  t >_  tg  > 0. 
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SIGNIFICANCE  AND  EXPLANATION 

\ 

In  this  paper  we  discuss  the  numerical  solution  of  a problem  which  arises  in 
polymer  rheology:  A mathematical  model  was  derived  to  descri're  the  elastic 

recovery  of  molten  plastics.  We  replace  this  model  by  a discrete  model,  for  which 
the  solution  can  immediately  be  constructed  by  computer. 

The  main  part  of  this  paper  concerns  the  relation  between  the  solutions  to  the 
original  and  discretized  models.  The  original  problem  is  a Volterra  equation; 

Volterra  equations  usually  occur  when  one  models  evolutions  which  depend  on  their 
history.  We  would  like  to  emphasize  that  when  one  solves  Volterra  equations  numeri- 
cally, it  is  important  that  the  discrete  problem  have  similar  qualitative  properties 
to  the  original  one.  In  the  particular  problem  considered  here  this  forces  us  to 
choose  a rather  poorly  convergent  method  if  we  wish  to  guarantee  the  results  obtained. 

The  numerical  results  confirm  a discrepancy  between  theory  and  experiments  indi- 

/ 

cated  previously  in  {3] , namely  that  when  the  elongation  of  a filament  is  large  and 


rapid,  the  model  predicts  somewhat  more  recovery  than  is  observed  experimentally. 
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The  responsibility  for  the  wordinq  and  views  expressed 
lies  with  MRC,  and  not  with  the  author  of  this  report. 


in  this  descriptive  summary 
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NUMERICAL  SOLUTION  OF  A SINGULARLY  PERTURBED 


NONLINEAR  VOLTERRA  EQUATION 

Olavi  Nevanlinna 

1.  Introduction. 

In  this  paper  we  consider  the  numerical  solution  of  some  singularly  perturbed  nonlinear 
Volterra  integrodifferential  equations.  These  were  modelled,  after  a problem  arising  in  poly- 
mer rheology,  and  studied,  by  Lodge,  McLeod  and  Nohel  [3], 

The  solutions  to  these  equations  have  a boundary  layer  for  small  positive  time  values 
but  one  is  actually  mainly  interested  in  the  behavior  as  time  is  large.  The  solutions  are  not 
exponentially  stable,  so  that  errors  during  the  time  stepping  generally  affect  the  estimated 
behavior  at  infinity.  However,  we  show  that  after  any  positive  time  * ^ , the  convergence 
for  the  scheme  we  propose,  is  of  order  O(h^)  with  some  S e (0,1),  uniformly  on  t >_  tQ 
and  0 < y <_  y 0 , where  u is  the  singular  perturbation  parameter.  On  any  compact 
[tQ,  T]  C (0,“)  the  convergence  order  is  0(h),  uniformly  again  for  u e (0,  uQl  . 

In  order  to  obtain  such  a strong  result  the  discretization  has  to  be  chosen  to  preserve 

as  much  properties  of  the  continuous  problem  as  possible.  This  we  achieve  using  the  backwards 

nh  n 

differences,  i.e.  we  replace  x’  (nh)  by  h 1(x(nh)  - x((n-l)h))  and  / rtslds  by  h£  <?(  jh)  . 
It  seems  unlikely  that  similar  results  could  be  proved  for  the  second  order  central  difference/ 
trapezoidal  rule  discretization.  In  fact,  in  order  to  prove  results  independently  of 
he  (0,  Uq)  we  have  to  preserve  the  pointwise  monotonicity  of  the  discrete  solutions,  other- 
wise the  nonlinearity  can  cause  them  to  run  out  of  control.  As  a very  simple  model  problem, 
consider  the  numerical  integration  of 

(1.1) 

using  one-step  methods 


hx'  + x = 0,  t > 0,  x(0)  = 1 


—(x  - x ,)  + 0x  + (l-8)x  , = 0 

h n n-1  n n-1 

These  methods  preserve  boundedness  uniformly  in  (J  > 0 iff  9 2.  1/2,  and  pointwise  monotonicity 
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when  0 > 1.  Since  for  0 _>  1/2  the  local  truncation  error  increases  with  0 it  is  reason- 
able to  have  0 = 1,  i.e.  the  backwards  difference  method. 

For  ordinary  differential  equations  it  is  possible  [6]  to  derive  sharp  error  bounds  if 
we  know  for  which  Ah  e (t  the  method  (linear  multistep  method)  yields  bounded  solutions  to  the 
test  equation 

(1.2)  x*  = Ax  . 

The  set  of  such  values  of  Ah  are  customarily  said  to  form  the  stability  region  of  the  method. 
A similar  concept  for  Volterra  integrodif ferential  equations  can  be  based  on  a test  equation 
of  the  form 

t 

(1.3)  x'  = ax  + 8^/  x(s)ds  , 

(or  similarly  for  integral  equations  of  the  second  kind),  see  (1),  [2).  However,  this  test 
equation  is  essentially  an  autonomous  differential  equation,  and  one  cannot  hope  that  preserv- 
ing boundedness  of  solutions  for  such  a test  equation  would  imply  good  error  behavior  for  pro- 
blems where  the  kernel  is  truly  time  dependent.  In  [4)  we  proposed  a different  type  of  test 
for  discretizing  the  integral  term:  often  in  applications  the  kernel  k(t,s)  satisfies  some 
positivity  properties  of  the  form: 

T t 

(1.4)  / >£(t)  / k(t,s)  v>(s)ds  dt  ^ 0 

0 0 

for  all  T > 0 and  all  smooth  <f.  We  characterized  those  discretization  schemes  which  pre- 
serve this  kind  of  positivity,  and,  in  [5),  we  derived  error  bounds  for  such  methods  when 
applied  to  some  systems  of  nonlinear  Volterra  equati  s.  We  also  applied  the  technique  to 
ordinary  differential  equations  and  showed  that  for  linear  multistep  methods  preserving  the 
positivity  was  equivalent  to  A-stability  (the  left  half-plane  belonging  to  the  stability 
region) . 

For  the  problem  studied  in  this  paper,  the  global  stability  of  the  solution  is  not  a con- 
sequence of  a positivity  property  of  the  form  (1.4),  but  it,  instead,  follows  from  a pointwise 
positivity  of  the  kernel  and  pointwise  monotonicity  properties  of  the  nonlinear  term  and  the 
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solution  itself.  Therefore,  we  asked  that  these  should  be  preserved  in  the  discretization. 
That  narrowed  the  class  of  method  down  to  the  first  order  backwards  difference  method,  but  in 
return  we  can  prove  a fairly  strong  convergence  result. 

We  state  our  theorem  in  section  2,  prove  it  in  section  3 and  discuss  the  application  to 
polymer  rheology  in  section  4. 

Finally,  1 would  like  to  thank  Professor  J.  Nohel  for  introducing  me  into  this  problem. 


2.  Convergence  theorem. 

The  following  integral  equation  was  studied  in  (31: 


(2.1) 


U*'  (t)  + / a(t-s)F(x(t),x(s))ds  = 0,  t 0 , 


under  the  assumptions: 


(H  ) 
a 


x(t)  = g(t) , t < 0 


a c 10, <»)  Ci  L^tO,*),  aft)  > 0,  a'(t)  < 0 


^a’(t)/a(t)  nondecreasing  , 

C F : R+  x ]R+  ■*  R,  F(x,x)  =0  for  all  x > 0 , 


( H ) 


'A'  ^ 


I 

(H  ) < 


C (R+  X K+)  and  Fj(y,z)  > 0,  F.,(y,z)  < 0 


(y,z  • R+  ; subscripts  denote  partial  differentiation) 


,'g  : ( — ,0)  - (0,»);  g{— ) = 1 , 


g |^q  is  nondecreasing  and  g(0)  > 1 

We  state  a convergence  theorem  for  the  following  discretization  of  (2.1): 

n 

“(x  - x ) + h } a(nh-jh)F(x  ,x.)  = 0,  n > 0 

h n n-1  J n l 

1 = — od  J 


(2.2) 


3 = 


x_.  = g(  jh)  , j 0 . 


We  shall  assume  that  the  initial  function  q is  also  differentiable: 

^ q satisfies  (H  ) and  q > C*(  — ,0]  , 
i 9 

(H')( 

q sup  | g • (s) I < « . 

Sf  IP 


Theorem.  Assume  that  (H  ),  (H_)  and  ( H * ) hold.  Then  for  all  u > 0.  h ' 0.  conation 

<1  F cj 

(2.2)  has  a unique  solution  (x  1,  such  that  x . (l,q(0)l  and  for  n > 0,  x < x 

n n n - n-1 

There  exist  constants  (j0  ■ 0,  hn  ■ 0,  Cj,  C2>  C^,  C4  > 0,  C^,  such  that  for  any  T < «■, 


4 


/■  yx'  (t)  + y>(x(t) ) - J a(t-s)  ^>(x(s) ) ds  = f(t),  t _>  0 , 

l 

v x(0)  = g(0)  , 

0 

where  f(t)  = J a(t-s)  ip(g(s) ) ds.  In  particular,  from  this  equation  we  see  immediately  that 
the  limit  equation  (y  = 0)  is  numerically  well  conditioned,  as  implied  by  (2.3). 


-5- 
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3 . Proof  of  the  theorem. 


In  order  to  prove  our  theorem  we  have  to  first  show  that  the  discrete  equation  has  essen- 
tially all  the  properties  that  the  continuous  problem  has  and  in  particular  all  bounds  hold 
with  constants  independently  of  the  step  length  h. 

Denote  a(nh)  = a , then  {a  } e a >0,  Va  = a - a < 0 and  (~)  is  nonde- 

n n n n n n-i  n 

creasing.  We  consider  the  discrete  Volterra  equation 

n 


(3.1) 


Vx  + hV  a . F(x  ,x.)  = 0,  n > 0 
h n L n-j  n ] 


x . 
3 


for  j < 0 


Here  g_.  = g(jh)  so  that  lim  g^  = 1 

Lemma  1.  There  exists  a unique  {x  } 
n > 0. 


and  { g . } 

J 

satisfying 


is  nondecreasing  and  gQ  > 1. 

(3.1).  Moreover,  1 < x < g. 

n u 


for  all 


□ 


Proof.  Write  (3.1)  into  the  form 


(3.2) 


l 


a 


n-3 


F(x  , x . ) 
n 3 


x 


1 


The  left  hand  side  of  (3.2)  is  a continuous  increasing  function  in  x^,  mapping  [ 1 » onto 
[a,b] , where  a < 1,  b > gQ.  Therefore  it  has  a unique  solution  x^,  in  the  open  interval 
(l,g0)  , whenever  x^  e [l#g0l. 

□ 


Lemma  2.  {x  ; is  nonincreasing  for  n > 0. 
n — 


Proof . 

By  Lemma  1 

x,  < x . Let 

1 0 

m > 2 

be  first 

index  such  that  Vx  >0. 

m 

We  show  that  then 

2 

V x < 

m — 

0,  being  a 

contradiction . 

Since 

F(x,x)  = 

0 , we  have 

rn- 1 

£ V2x  + 
h m 

h ) a 

^ m- 

. JF ( x ,x . 
-3  m j 

) - F(x  ,x . ) ] 

m—  i j 

(3.3) 

m-2 

= -h  ) Va  . F(x  . ,x . ) 
L m~3  m-1  i 


-6- 


By  assumption  so  that  the  second  term  on  left  in  (3.3)  is  positive.  If  we  can 

show  that  the  term  in  the  right  is  nonpositive  then  V^x  must  be  negative  and  a contradiction 

m 

follows . 

Define  J to  be  the  largest  index  j such  that  g.  < x . Set 

j — m— 1 

m-2 


and 


* h Y Va  .F(x  ,x . ) 
L m- j m-1  i 

— OO  J J 


R * h y a . F ( x , X . 
m L m-j  m-1  j 


Then  R is  negative  which  follows  from  the  identity 


m-1 


Rm  = * h Vxm  -hJa  {F(x,x.)-F(x  , ,x. 

m " m m-3  m 3 m-1  3 


)} 


and  from  the  assumption  that  x > x . But 

m m-1 


Va 

s„  > s - (— ) R 

m — m a m-J  m 


= . - <22* 

a m-3  a m-J  j 

m-2 

+ h l { (— ) . - <^£)  }n.  , 

j+!  a m-3  a m-J  3 

where  rn  = am_jF(xm_1,x;. ) . Since  (^)  is  nondecreasing  the  term  in  brackets  in  the  first 
sum  is  nonnegative  and  in  the  second  sum  it  is  nonoositive.  But  by  construction  n.  is  non- 
negative for  j < J and  negative  for  j > J and  therefore  S >0,  and  a contradiction 

m — 

follows . 

□ 

Next  we  formulate  some  other  properties  of  the  solution  {x  },  all  of  which  have  counter- 

n 

parts  for  the  continuous  problem,  [31. 

Proposition  1.  Solution  depends  monotonically  on  the  initial  data:  if  q < q for  all 

3 - j 

j < 0,  then  x < x for  all  n > 0. 

— n — n 

Proof . Let  m be  the  first  index  such  that  x > x . Then 

n n 

-7- 
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-'Vx  - Vx  ) = 
m m 


-h  ) a . ( F(x  , x . ) -Fix  ,x . ) 
_ m-j  m j m j 


+ Fix  ,x.)  - Fix  ,x.)  } 
m j m j 


and  since  Fix  ,x.)  > Fix  ,x.)  and  F(x  ,x.)  - Fix  ,x.)  we  get  Vx  < 7x  which  is 

rag—  mj  raj—  mj  m — m 


tradiction. 


Proposition  2.  Solution  depends  monotonically  on  y:  if  {x  (y  )),  {x  (y  ) ! are  two  solu- 

n 1 n 2 

tions  to  (3.1)  with  the  same  initial  data  but  different  y’s,  then  y^  - y^  implies 

x I m , ) < x (y  ) for  all  n > 0. 
n 1 n 2 


Proof.  One  sees  easily  that  x^(y^)  < SuPPOSG  therefore  that  m is  a first  index 

such  that  x (y  ) > x (y  ) . This  would  imply  V(x  (y_)  - x (y,))  <0.  But 
m J-  — m 2 m2  ml 

h m2  ml 


~~  I am  .{F(x  (y  ),  x . (y  ))  - Fix  (y  ) , x.ly  ))} 
y~  m-j  m2  J 2 ml  -j  1 


+ [■ — ] h l a . F ( x (y  ) , x.(y  )) 

U i V2  m-3  ml  j 1 


which  shows  that  the  right  hand  side  of  (3.4)  is  nonnegative.  In  fact. 


=-TTVxm(V  -° 


Fix! y ) ,x . (y  ) ) - F ( x ( y ) , x , ( y ) ) 
m 2 j 2 raljl 


= {Fix  (y  ),x.(y  ))  - Fix  (y  ),x.(y  ))) 
m 2 32  m 1 32 

+ {Fix  (y,),x.(y  ))  - Fix  (y  ),x.(y  ))}  < 0 , 

ralj2  mljl  — 

since  both  bracket  terms  are  separately  nonpositive.  Hence  we  arrived  into  a contradiction. 


4 


n-i 


Then 


0 < S - h l (— ) .a  . F ( x , X 

, a n-i  n-i  n-1 


n-1 


- (T>ih  I VjF(Vi-  V lcl(X>ilh  l aj 

1 n 


|(~>iIy{/  a(s)ds  + OlhlJjVx^l  , 


Va 


since  — is  negative  and  nondecreasing. 


F(x  , x.)  < 0 , and 

n-1  D - 


n-1 


-h  I a 1F<X„  , < x. ) = rVx  + hi  a ,F(x  ,g .) 
“ n-]  n-1  ] h n L n-g  n 

1 —CO  J J 

n-1 

-h  l a .F  ( 5 . ,x . ) Vx 
1 "'I  1 1 D n 


so  that  C = sup  F(y,z) 
y,ze [l,g0l 


For  the  other  sura  we  get 


,1  < h£-(— ) .a  . | F(x  , g.) 
1 - „ a n-g  n-]1  n-1'  yj' 


< c|(^>  |h  J1  a. 
an1  L 3 


and  together  they  imply  (3.7). 


n ,n  - - h yn  + Gnyn  = h Yn-1  + fn 


Set  Vxn  = yr>*  T^en  we  obtain  from  (3.6)  that  {j- 
and  further 

(h  + Gn) lvnl  -hlyn-ll  + lfJ  • 

Suppose  from  here  on  that  h <_  h^,  where  h^  is  so  chosen  that  1 - | (— ) | > 1/2 

we  see  from  the  bounds  for  G and  f that  for  some  constant  G > 0 


Then 


n n 


(3.8) 


(h  + G)lyJ  - hlyn-J 
-10- 
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where  g < c|(— ),|h  Y a..  Since  (— ) , = * h + Olh^)  we  obtain  from  (3.8) 

n — ' a 1 1 ^ j al  a(0) 


lyj 

H i=2  k=j 

u J 


Consider  the  second  term  in  the  right  hand  side  of  (3.9): 


S = l (1  ♦ — ) h l a, 
n \i  , . k 


By  partial  summation  we  obtain 


2]  n-l(l  + — > -1  (1  + — ) -1  ~ > 

1 + M + hnTTi h I • 


U r . "2  hG/u  j hG/u 


But  a . = 


'j  ” " ? - ’(  a'i  l \ 

k=g+l  k=3+l 


2 n-1  (1  + — ) -1  m ( , Va  . 

,,  . hG,  v U a . < — ( — ) , S 

( 1 + — ) ) — — — 1 — hG ' a 1 ' n 

V >2  hG/u 

and  hence,  if  u , and  h <_  h.  are  such  that  1 - £ ^-|  (— ) | > 1/2  for  u € (0,u_] 
u uj.  G h a i — u 

7a  a' (0)  2 

h € (0,h  1 (which  is  possible  since  ( — ) = — — — h ♦ 0(h  ) ),  we  finally  have 
0 a 1 a ( 0) 


._  n+l  . 2 «• 

S < 2 [(1  + ^)  - (1  + — ) ]h  l a . 

n = hG  U u *•  k 

k=n 


|y:l  - \\  - 90l  ic  J h l a.  , 


(3.5)  follows  by  substituting  the  bound  (3.10)  for  S into  (3.9). 


The  next  result  will  be  used  to  Drove  the  uniform  convergence  of  the  approximations  on 


I 


the  half  line  t > 0. 


Corollary  1.  If  additionally  { ja^ } e 5,  , then 


< C (1  + — ) + C h l (k-n)h 

1 w k=n+l 


Proof . (3.11)  follows  immediately  from  (3.5). 


Up  to  this  point  we  have  more  or  less  followed  the  treatment  of  the  continuous  problem  in  [3) . 
In  order  to  prove  convergence  we  need  to  investigate  the  perturbation  properties  of  the  differ- 
ence equation.  When  proving  the  convergence  we  would  like  to  substitute  the  local  truncation 
error  in  the  place  of  the  perturbation. 

However,  we  can  control  the  solutions  only  when  we  allow  nonnegative  perturbations,  and 
the  local  truncation  error  does  generally  change  signs  during  time  stepping,  so  that  some  extra 
consideration  is  in  order. 

Lemma  3.  Assume  u.  = v.  = g.,  j < 0 and  {u  },  {v  } satisfy 
3 3 J - n n 


q-  Vu  + h y a ,F  (u  ,u.)  = a 
h n L n-3  n j n 


(3.13)  £ Vv  + h l a F(v  ,v  ) > q . 

h n oon“^  n ] — n 

Then  v > u for  all  n > 0. 
n — n — 

Proof.  Let  m be  the  first  index  such  that  v < u . Then  V(v-u)  < 0 and  we  get 
mm  m 


It  V(v-u)  + h y a ,{F(v  , v . ) - F(U  , u . ) } < 0 , 

h m L m-]  m 3 m j 


contradicting  (3.12),  (3.13). 


Let  now  {y^}  be  a sequence  from  which  it  is  known  that  fnr  all  n yn  t (1,0^1  (we  know 


that  x(nh)  c ( 1 , <3q1  ) an<5  that  it  satisfies 


-12- 


r 


11 

r-  Vy  + h F a .F(y  ,y.)  = p , n > 0 
h n L n-j  n j n 


and  y.  = g.,  for  j < 0.  Construct  two  seauences  {w  } and  (z  } as  follows:  w satisfies 
]]  ~ n n n 


j"  Vw  + h J a .F(w,w.)  = n 
h n n-i  n y n 

where  w.  = g.  for  j < 0 and  a = |p  [ if  this  will  imply  w *■  g . otherwise  set 
13  nn  n—0 

w = g and  compute  the  corresponding  value  for  q . The  other  sequence  (z  t will  serve 
n u n n 

as  a lower  bound.  Therefore  define  it  by  the  same  equation,  now  with  perturbation  r^  = -|d  ! 
whenever  you  would  have  >_  1,  otherwise  set  = 1 and  define  correspondingly. 

Lemma  4.  If  {x  } satisfies  (3.1)  and  (v  },  {w  }.  { z } are  aiven  as  above,  then 
n n n n 

I x -y  I < I w -z  I . 

n n — n n 1 □ 

Proof.  By  construction,  q > 0 and  r < 0 so  that  bv  Lemma  3 we  have  z < x < w . The 
n — n—  n — n — n 

conclusion  will  follow  since  we  also  have  z < y < w . To  see  this,  assume  for  example  that 

n — n — n 

m would  be  the  first  index  for  which  y > w . Then,  since  y < we  must  have  q = Id  I , 

mm  ' m — ^0  n 1 ' n 1 

and  therefore  we  are  back  in  the  situation  of  Lemma  3 and  a contradiction  would  follow. 


According  to  Lemma  4 we  can  proceed  as  follows.  Suppose  that  x(t)  is  the  solution  to 
the  continuous  problem.  Define  the  local  truncation  error  sequence  {t*1}  by 


n 

^■(x(nh)  - x((n-l))  + h ] a .F(x(nh),  x(jh))  *=  t 
h n-3  n 


If  x satisfies  (3.1)  with  g.  * g(jh),  then  in  order  to  bound  lx  - x(nh) I it  is  suffi- 

n 3 1 n 1 

cient  to  bound  |w  - z I,  where  {w  - z } satisfies 
n n n n 


ii  , 

0 < ^-V(w-z)  h 7 a .{F(w  , w . ) - F(z  ,z.)}  < 2 1 t | 
- h n L n-3  n j n j = 1 n 1 


with  Wj  * Zj  = g^ , j <_  0,  and  by  the  above  construction  we  can  assume  that 

1 < z < w < for  all  n > 1. 

— n - n—0  — 
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Lenmd  5.  Assume  that  (w  I,  (z  } C [l,q  ] and  satisfy 
n n 0 


(3.16) 


0 < ^-V(w-z)  + h y a . (f(w  ,w.)  - F(z  ,z.)l  < p , 

— h n ^ n-j  n j n j = rn 


and  w.  = z.  - g^  for  j <_  0.  If  yQ,  hQ  are  small  enouqh,  y € (0,uQ) , and  h f (0,hQl , 
then  for  some  constants  A,  B,  C,  depending  on  yQ  and  h^  but  not  on  y and  h we  have 


(3.17) 


Bh 

v 

|w  -z  | < Ch  T p e~  * -^,U-  "j1  p (1  + ^)‘J 
n n “ v=0  n’V  y=0  n_V  W 

y 


Proof.  As  in  the  proof  of  Lemma  3 we  observe  that  (3.16)  implies  w > z for  all  n ■>  0. 

n — n — 

Therefore  we  can  devide  the  sum  in  (3.16)  into  two  parts: 


n-1 


(3.18) 


r~  V (w-z)  + h T a ,(f(w  ,w.)  - F(z  ,w.)) 

h n J-  n~]  n ] n ] 


n-1 


P + h J a ,{F(z  ,z.)  - F ( z ,w.)} 
n v n-]  n i n i 


Let  a = min  F (£,n)  and  8 = max  -F„(£,n)  , then  (3.18)  and  w.  > z imply 


C.ni  [l,q 


(3.19) 


[ 1 , 9q] 


n-1 

Jc  V (w-z)  + ah  y 

n n “ 


a . • (w  -z  ) 
n-j  n n 


n-1 


< p + Bh  y a . (w.-z.) 

- n j n-j  i i 

OO  00 

Set  w-z  = C and  Ba(0)  = B.  Choose  A < aj  a(s)ds  and  h so  small  that  ah  £ a(jh)  > A. 

0 1 

Then  (3.19)  implies 


(3.20) 


where 


(i+— k <t  + — h y z,  + 

y n = n-1  y J j 


n-1 


— p , n > 1 
Vi  n 


Ah  ,_1  h 


C,  < (1  + — ) - p, 
1 - y y *1 
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To  solve  (3.20)  we  first  consider  the  roots  of  the  polynomial 

(1  + — >\2  - (2  + ^ + 1 . 
U p p 


When  Ug  and  hg  are  small  then  we  have  after  a straightforward  estimation  two  positive 
roots,  Aj,  *2  say,  such  that  for  p * (0,poJ,  h ( <0,hg] 

1 + - -Kh2 

'^-rrsr-1 


and 


l2i' 


Bh 

A 


where  K depends  on  p and  h„ . The  solution  C then  satisfies 

u u n 


(3.21) 


hi  1 n'1 

n ~ U a .Ah  A -A  ^ _ ^n-v^ 


1 ♦ — 2 1 v=0 

u 


NOW  V‘iirft  and  ui-m2ir(lt7,eA'  and  -(Xi  - + “>”■ 

(after  possibly  choosing  a smaller  A and  hg) . Substituting  these  bounds  into  (3.21)  yields 


V+1  _ Av  . Av+1  + 

2 2 Ai  *r 


Bh 


n-1 


- , 1 V , Bh  . Bh. 

Cn  - A Pn-v(T  (1  + T> 


Bh 

Bh.  A Ah/p 


v=0 


n-1 


< i + i 

l+Ah/p  1 p ’ ’ 


Bh 


Ch  l e'p 


v-0 


n-v 


, . n-1  -v 

Vm  c M , Ah 

l+Ah/p  p ' pn-v 


for  some  constant  C,  and  we  hence  obtained  (3.17) 


Next  Lemma  bounds  the  local  truncation  error. 
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Lemma  6.  Let  { denote  the  local  truncation  error  (see  (3.14)).  Then  there  exists 

c = c<y  ,h  ,q  ) and  K = K*uo'qo*  such  that  for  U f (0#W0J  # h * (0'V  We  have 

I h,  ^ _ ( r . ^ 1 -K(n-l)  h/W 

t < Ch(  J a s ds  ♦ - e 

1 n 1 - ' i , U 

\ nh-h 

(3.22) 

+ / a(nh-s)£/  a(s)ds  + - e , 

0 s U ) 

□ 

Proof.  If  x(t)  is  the  solution  to  the  continuous  problem,  then 

n jh 

(3.23)  lix’(nh)  + £ / a (nh-s)  F ( x !nh)  , x(  s) ) ds  = 0 

(j-l)h 

Let  us  consider  the  sum  first: 

n jh 

£ / a(nh-s) F(x(nh) ,x(s) ) ds 

-°°  (j-l)h 


n 

= h J a(nh-jh) F(x(nh) ,x ( jh) ) + 

where 

n jh 

S = 1 / ( ( j-1) h-s) [-a' (nh-s) F(x(nh) ,x (s) ) ds 

-»  (j-1) h 

and 

n jh 

S = J J [ ( j-1) h-s] [a(nh-s) F (x(nh) ,x(s) ) ] x' (s) ds 

-»  (j-l)h 

We  know  that  x(t)  is  decreasinq  and  x(t)  t (1,9^1  for  t > 0 [Theorem  1:3].  Therefore, 

set  C = sup  |F(y,z) ! and  D = sup  ]f  (y,z) | . We  split  S. 's  into  two  parts  and 
y,zc[l,g0l  y,z<-[l,qnl 

estimate : 

0 jh 

| S | = |£  J [ ( j-n)h-s] [-a’ (nh-s)F(x(nh) ,x(s) ) ]ds 
1 — (j-l)h 
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0 jh 

<_  Ch  £ / |a'(nh-s)|ds 

(j-l)h 


-c|Tw"|h^  a(s)ds  1 


nh 


where  we  again  used  the  logarithmic  convexity  of  a(s).  Similarly  we  get,  since  x(t)  is 
decreasing, 

, . . nh 

|S,, | <_  h|  "■  ■ --|  / a(nh-s)  (-F(x(nh,x(s) ) ) ds 
nh  0 

But  -/  a(nh-s)F(x(nh)  ,x(s)  )ds  = ux'(nh)  + / a(nh-s)  F (x (nh)  ,g ( s)  ) ds , and  since  x'(nh)  *■  0, 
o -® 

we  obtain 


lsnl  ^clT(f-lh  J a(s,ds  * 


nh 


Hence 


lsil  - 2clT<mLlh  ! a(s,ds  • 

nh 

Also  is  of  this  form: 

0 jh 

|s  | =|£  j [ ( j-1) h-s) a(nh-s)F  (x(nh) ,g(s) )g ' (s) ds I 
-»  (j-l)h 

00 

£D  sup  g’(s)-h  j a(s)ds  , 
se(-°°,0]  nh 

while  S^2  carries  the  effect  of  the  boundary  layer:  By  [Theorem  2:31  (which  should  be  com- 
pared to  Proposition  3)  there  exists  constants  and  K0  such  that 


(3.24) 


K -K  t/U 

0 < -x' (t)  < — e + K.  f a(s)ds 

- - w 2 i 


Substituting  this  into  S^^  yields 


nh 


k2  -KjS/u 


22 


| < Dh  / a(nh-s)  ( — e + K j a(r)dr)ds 
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We  have  to  still  consider  the  error  introduced  by  replacing  Px*  by  ^ Vx.  By  mean  value 


theorem  we  have  for  some  F,  ((n-l)h,nh) 

n 


( 3 . 25) 


x(nh)  - x((n-l)h)  . , . . 1 , 4 

p : r nx'  nh  - t hx'  U ) 

h 2 n 


Using  line  (4.1)  in  (31  and  (3.24)  we  obtain 


' x"C  ) Ch  ■ 
2 n . 


-W*' 


♦ / a(s)ds> 

f ( 


Collecting  all  terms  yields  (3.22). 

n 

Note  that  (3.22)  implies  that  |t^|  is  uniformly  bounded  for  n _>  2,  g e (0#pgj , h a (0,hg] . 

However,  the  same  is  true  also  for  |t^'|  , which  we  now  demonstrate.  We  see  from  the  proof 

above  that  the  unbounded  term  in  the  bound  for  t!?:  C~,  comes  from  (3.25).  However,  the 

1 P 

solution  has  a boundary  layer  for  small  t 0,  and  this  implies  the  uniform  boundedness  of 

this  term.  In  fact,  there  exists  a continuously  differentiable  function  £ on  [0,®),  such 

that  f,(t)  • (l,g  | and  sup  ;f,'(t)l  < ®,  and  a constant  K = K(h  > , such  that 

tr [0,®) 


(3.26) 


'x(h)  - C (h/vO  \ Kh,  for  h c (0,h 


0 ' 


(see  (Theorem  7;  3]  and  notice  that  since  £.(t)  is  decreasing,  so  is  |£'(t)  |).  But  (3.26) 
implies  that 


X<h)-qn 


yx’  (h) 


5<h/y)-qn 
I 9 

1 h/y 


+ | V (h/y)  ! + 2Ky 


which  is  bounded  as  soon  as  y j_  y(),  and  therefore  there  exists  M = M(y0,hQ)  such  that 
|TJ'  < M for  y ■ |,0,  h : hQ. 

- gt 

Next  we  simplify  (3.22)  somewhat  by  assuminq  that  alt)  <_  ne  for  some  a > 0,  8 i 0. 

Then  we  have  for  some  C • « and  •,  0,  with  y ■ 0 if  B > 0, 


Ch(e*'(n*Uh  4 i e-K<h-l)h/y)(  n > j 


and 


C. 


-1(3- 


4 


Substitute  these  into  (3.17) 


The  first  term  in  the  right  of  (3.17)  yields, 


(3.27) 


C1eBT/Ah,  if  Y > 0,  and 


(3.28) 


C1(l+T)eBT/Ah,  if  y = 0 


The  other  term  yields 


and. 


and 


h (1  + AH,'1  [(1  + Ah,'n+1|h|  + "f  |Th  |(1  + Ah,'V\ 

p p ^ p 'l1  “ 1 n-v 1 P J 


= C - (1  + — ) + S,  where  S = S + S,  , 

P P 12 


. I < - (1  + — ) -ch  J Ml  + — ) < | h 

1 1 — VI  P „ P — A 


2 1-n  1-n 

Is  | = (h,  2 

1 21  V a - 6 


where  a = 1 + — ana  8 = Taking,  if  necessary,  a smaller  A,  we  see  easily  that 


Kh/p 


|S2|  is  bounded  in  the  form 


1-n 


Is  I <ck|(1^‘  (l.^h)1'”} 

1 2 1 — U I P P J 


Summarizing  the  above  bounds  yields  (2.3).  To  complete  the  proof  assume  that  a(t)  decays 
exponentially  so  that  (3.27)  holds.  From  Corollary  1 we  see  that  for  m ^ n,  and  for  some 
Y > 0, 

t I -Ynh  , , Ah.  n 

n m1  — p 

and  for  the  continuous  problem  we  have  similarly  (see  [Corollary  2 . 1 ; 3 ] ) 

|x(nh)  - x(mh) | < Ce 


Assume  that  mh  > T and  nh  < T. 
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Then 


x(mh)  - x < x(mh)  - x(nh)  + x(nh)  - x+  x - x , 
m — 1 11  n n m 

where  the  right  hand  side  is  bounded  by  terms  of  the  form 

C T 

C1he  2 , C3e'YT 

(and  by  the  boundary  layer  terms) . Allowing  T to  grow  with  h one  can  bound  these  both 
6 c -1 

terms  by  Kh  , with  6 = (1  + 2/y)  , which  completes  the  proof  of  the  theorem. 


4.  Example  from  polymer  rheology. 


For  background  material  see  (3). 

In  this  example 

x(t)  denotes  ■ 

plastics  filament:  it  is  originally 

of  length  1 

/ then 

it  is  elongated 

is  allowed  to  recover  freely.  Hence 

x(t)  = g ( t ) 

, t < 0 

describes  the 

x(t)  for  t > 0 the  free  recovery. 

Functions  a(t) , F(y,z)  and  g(t)  are  as 

follows 

8 

-t/l 

a (t)  = ) a. 

1 

e 

, 

i=l 

where 

i 

T . 

1 

a . 
i 

1 

io3 

10" 

3 

2 

io2 

1.8 

3 

10 

1.89 

• 10+2 

4 

1 

9.8 

• io3 

5 

io-1 

2.67 

• io5 

6 

io-2 

5.86 

• 106 

7 

io-3 

9.48 

• io7 

8 

io'4 

1.29 

• io9 

(these  numbers  were  provided  by  Professor  Lodge  and  based  on  results  of  H.  M.  Laun) , 


F(y,z) 


and 


9(t)  =< 


'L.  . « 


, t < -t„ 


K(t+t0> 

e > -to  i t i 0 


Since  a,  F,  and  g are  all  of  relatively  simple  form,  one  can  write  the  inteqral 

0 

/ a ( t-s) F (x , g(s))ds 
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down  exactly  and  only  discretize  the  other  oart 


J a(t-s)F(x,  x(s))ds 
0 

Furthermore,  although  the  backwards  differences  lead  to  an  implicit  method,  the  specific  form 
of  F implies  that  at  each  step  one  has  to  find  a real  root  to  a cubic  polynomial,  and  that 
can  be  written  down  in  terms  of  square  and  cubic  roots.  Therefore  the  final  recursion  formula 
is  essentially  an  explicit  method,  which  generates  from  the  finite  information  (x^,...,xn) 
a new  value  x 

n+1 

In  an  earlier  work  (see  (Appendix  A; 3])  equation  (2.1)  was  considered  with  y = 0,  and  it 
was  found  that  the  model  predicted  more  recovery  than  was  observed  experimentally.  Since  the 
solution  is  increasing  with  y,  introducing  y > 0 leads  to  a reduction  in  the  predicted  re- 
covery. Unfortunately,  the  numerical  experiments  show  that  the  final  recovery  is  effected  very 
little  by  y (for  y e [10  , 1]) , although  the  solutions  behave  very  differently  for  small 

t > 0,  and  so,  for  high  strains  there  still  is  a disagreement  between  theory  and  data. 
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